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Abstract: The superconducting transition in presence of strong columnar disorder parallel 
to the magnetic field is considered. A solvable model appropriate for description of the 
broad crossover regime towards the true "glassy" critical behavior is constructed, and the 
behavior of the thermodynamical quantities and of the Edwards- Anderson order parameter 
near transition is obtained. The critical exponents for the correlation lengths orthogonal and 
parallel to the magnetic field are the same and in agreement with the experimental values. 
The dynamical critical exponent is z = 2, also in agreement with the measured value. Several 
perturbations to the solvable model are considered and shown to be irrelevant for the critical 
behavior. It is argued that there exists an optimal density of defects at which the tansition 
temperature at given magnetic field reaches its maximum. 



1 Introduction 



The problem of superconducting transition in the magnetic field in presence of strong disorder 
has been a subject of great interest in recent years. From the point of view of applications, 
there is a need to know how to introduce defects into a superconducting material in a way 
that would maximize pinning of vortices and therefore increase critical currents On the 
theoretical side, understanding of the phase transitions in presence of quenched disorder and 
description of the low-temperature disordered phases have always been challenging problems, 
and often required novel physical concepts and mathematical techniques. A prime example 
of this is the theory of spin glasses. For the superconducting transition in the magnetic 
field, a variety of novel low-temperature phases have been proposed, differing in the cases of 
point like p|, |^ and line-like disorder [|, ||, |^ . In the present paper we discuss a theory of 
superconducting transition at high magnetic fields in presence of columnar (line-like) defects 
0. We have in mind a three-dimensional, anisotropic, strongly type-II superconductor 
(YBCO, for example) in typical magnetic fields of ~ IT, irradiated by a fiux of some heavy 
ions with energies ~ IGeV and with trajectories parallel to the external magnetic field. 
If the thickness of the sample in the direction of the beam is ~ 10/im, the ions are able 
to penetrate through the entire material, leaving behind the continuous tracks of damaged 
superconductor of diameter d ~ 50A. In the absence of disorder, the high-field fiuctuations 
of the order parameter, ^/'(r), are strongly enhanced by formation of Landau levels (LLs) for 
Cooper pairs. Such fluctuations lead to D ^ D — 2-dimensional reduction in the pairing- 
susceptibility, Xsci^jf), and, strictly speaking, eliminate the superconducting (Abrikosov) 
transition for D = 2, 3 0. The Abrikosov phase is then replaced by a new fluctuation-induced 
state, the density-wave of Cooper pairs (SCDW), in which the thermal average (iV'lr)^) 
a weak modulation not necessarily accompanied by a long-range phase coherence [0]. The 
chief effect of disorder is to remove the LL degeneracy and so to restore a possibility of a true 
superconducting transition. In this sense the superconducting transition that we will consider 
is induced by the presence of disorder. Xsc can now diverge at some flnite temperature, 
Tsc{H), determined by the strength of disorder. The superconducting transition corresponds 
to the Bose condensation of Cooper pairs into the lowest energy eigenstate of the random 
potential, which, we will argue, extends over the whole sample in situations of experimental 
interest. Furthermore, for experimentally relevant parameters, Tsc{H) can be far above the 
SCDW transition line P| over much of the H — T phase diagram, allowing us to treat the 
fluctuations that produce SCDW in an approximate way. 

In the following sections we flrst deflne a model for superconductor in a magnetic fleld 
parallel to the columnar defects, which represents a realistic description of the problem for 
certain ranges of magnetic flelds and densities of columnar defects. It is demonstrated that 
the model is exactly solvable. The solution exhibits "dimensional transmutation", i.e. the 
effective dimensionality of the transition changes continuously as a function of magnetic 
fleld. This effect is a direct consequence of analytic properties of LL wavefunctions and is a 
signature of the high- fleld limit. We determine the transition line m H — T phase diagram, 
the Edwards-Anderson order parameter, and the behavior of speciflc heat and magnetic 
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susceptibility in vicinity of the transition. There are similarities between the transition 
considered here and the one in the spherical model for spin-glasses 0. Using the information 
on diffusion in strong magnetic field that comes from the studies of the quantum Hall effect, 
we determine the critical exponents for the correlation lengths and the dynamical critical 
exponent. The calculated exponents are in good agreement with the experimental values 
01 . We then consider perturbations to our exactly solvable model. In particular, we discuss 
the breakdown of the model for higher density of defects. It is argued that disorder then 
becomes less efficient in removing LL degeneracy, enabling fluctuations again to suppress 
superconducting transition temperature, similarly to what happens in the homogeneous case. 
This naturally suggests the existence of the optimal dose of irradiation which produces the 
highest transition temperature at given magnetic field. We conclude by reviewing the salient 
points of our model and by briefly reanalyzing the approximations built into it. 



2 Statement of the problem and the solvable model 

We are interested in strongly anisotropic layered superconductors described by the Ginzburg- 
Landau (GL) Lawrence-Doniach model, with magnetic field perpendicular to the layers. 
Fluctuations in magnetic field are neglected (GL parameter k 1). We focus on the high- 
field limit, where the LL structure of Cooper pairs dominates the fluctuation spectrum: This 
is the case for fields above ~ {9/16)Hc2{0)(T/Tc{0)), where 9 is the Ginzburg fluctuation 



parameter |T0|. (For instance, in BSCCO 2:2:1:2, 6 ^ 0.045 and ^ 1 Tesla.) In this 
regime, the essential features of the physics are captured by retaining only the lowest Landau 
level (LLL) modes. Generally, the partition function is Z = J D[ip*,ip] exp(— S"), and 

b ^ 



^ = ^ E / d'f^vlM^ - V^n+i(r)r + 

n=l 

{a'{T,H) + XY.V{\f-r-imnm' + ^\Mm. (1) 
i ^ 

where a'{T,H) = a{T — Tc2{H)), b is the effective layer separation, n is the layer index 
and a, (3 and rj are phenomenological parameters. The magnetic field is assumed to be 
parallel to columnar defects, which are modeled by an effective potential V{r — fi) > 0, 
peaked at r = with a width comparable to the diameter of the columns d. The parameter 
A > represents the effective strength of disorder, and its positivness reflects the fact that 
the superconducting temperature is locally suppressed by damaging the material. Random 
variables in the problem are two-dimensional coordinates of defects, {i^}. We assume that 
the positions of columns of damaged superconductor are uncorrelated, i. e. that they are 
distributed according to the Poisson distribution Pn^tI, ...r^) = {e^P^p^)/N\ where P/v is 
the probability for finding impurities at the positions fl, ...r]v, A is the area of the system 



and p is the concentration of impurities [ITT 



The partition function in Eq. 1 is quite general, and defines the problem we want to 
study. To proceed, note that in certain range of parameters one may make two simplifying 
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assumptions. First, if the magnetic length / and the average distance between defects are 
both much larger than the diameter of the columns we may assume that V{r) = 6{r). At 
fields of ~ IT, / ^ 200A , and for moderate doses of 10^ — lO^^ion/cm'^ these conditions are 
reasonably satisfied. This approximation makes the random potential which describes the 
disorder completely uncorrelated in space, and facilitates the exact treatment. Second, after 
rescaling the fields and the lengths as {2bp27rP /TY^^iJj — ^ ^, r/(/-\/27r) r, the quartic 
term can be rewritten as 

d'r-\^Pn\' = ^EMn)i I rfV>„r)^ (2) 

where l3A{n) = {N / |'?/'n|^)/(/ is the generalized Abrikosov ratio corresponding to 
configuration ipni'f^, and N = A/2ttI'^ is the degeneracy of the LLL. We now observe that 
f^Ain) is only weakly dependent on the actual configuration, the well known example being 
the small difference in (3a between triangular and square lattice of zeroes [0. Thus, we 
may approximate f^Ain) in the quartic term by a constant of order unity. This replaces the 
local quartic term in the general theory (1) by an interaction still diagonal in layer indices, 
but infinitely ranged within a layer. This approximation neglects the lateral fluctuations 
that produce the SCDW transition [|^, and without disorder this theory has no finite tem- 
perature phase transition below four dimensions. Thus, the disorder is assumed strong on 
the fine energy scale set by the local quartic term, but weak compared to the energy scale 
of LL separation, so that the LLL approximation is sensible. The neglect of weak lateral 
correlations is justified if Tsc{H) is far above the SCDW transition line p. In that case 
the SCDW fluctuations start to matter only very close to the glassy transition and can be 
ignored in most realistic situations. Since superconducting and SCDW transitions arise from 
two distinct mechanisms, the respective transition lines scale differently in the H — T phase 
diagram and, for moderate disorder, we are assured of a wide crossover region near Hc2{T) 
where the neglect of SCDW fluctuations in our model is justified. 



3 Transition line and the order parameter 

After the /3a (^) — * {Pa) ~ 1 substitution the thermodynamics of the model becomes exactly 
solvable P]. We first introduce variables {x„} to decouple the quartic term and integrate 
over the fields This leads to Z = J H^dxriexp^—NS'), where 

5' = -E7/y+/ dVpf{V)Tr^^,m) 

n \Pa) jo 
n,m ^n,m—l ^n,m+l 

We drop the terms coming from the rescaling of ipni'f^) and introduce dimensionless combina- 
tions of GL parameters a = A/27r/^} x ^ (birP) / (T P) . The density of states for an 
uncorrelated random potential in the LLL can be found exactly by using the supersymmetric 
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formalism For the Poisson infinitely short-range scatterers it is given by the integral: 

Pf{V) = -Im^M r dte^^iiVt -f f ^(1 - e-^y)), (4) 
vr rfl/ JO y 

where / = p27rZ^ = Hij,/H, and is the matching field, at which there is precisely one defect 
per unit of flux. In the thermodynamic limit N ^ oo, the partition function is completely 
determined by the saddle-point of S'. Assuming that the saddle point is at x independent 
of the layer index, we finally write the free energy above the critical temperature 

F x"^ 1 

iViVzJ " ~ ^ 2 dVpf{V)\n[g,e{k)+ga + x + gxV], (5) 

where e{k) = 1 — cos(/c) and x is determined by the solution of 

4 J-i Jo gne{k) + ga + x + g^V 

In Eq. 6 it is important to know the behavior of density of states at low energies. For / < 
1, density of states has a delta-function singularity at \^ = 0, while for / > 1, PfiV) ~ V^'"^ 
when ^ |T3|. The transition line, Tsc{H), in the H — T diagram is determined by Eq. 



6 and x + ga = which corresponds to condensation of Cooper pairs into = and = 
eigenstate of the random potential. Without disorder the LLL is completely degenerate and 
piy) = S{V). The integral on the right hand side of Eq. 6 is then infrared divergent when 
X + ga = 0, and there can be no finite temperature phase transition. With the columnar 
disorder present, from the behavior of the density of states it is easily seen that there will 
be a non-zero transition temperature only if concentration of impurities and magnetic field 
are such that / > 3/2. Below this value of / LLL degeneracy is not sufficiently lifted by the 
random potential and thermal fluctuations prevent a finite temperature phase transition in 
our model. / = 3/2 determines the effective lower critical dimension for our model. After 
introducing dimensionless quantities t = T/Tc{0), h = H/Hc2{0) and A' = XHc2{0)/(f)oaTQ, 
where 0o is the fiux quantum, we perform the integration over wave-vector k in Eq. 6 to 
obtain the expression for transition temperature: 



tsc{h) = {l-h) 



2A' Jo Jv^ + {2r]V)/{hX'aTo 



(7) 



Notice that when A' we have tsc{h) 0, while for increasing A', transition temperature 
tsc{h) increases. Also, with increasing density of columnar defects (increasing parameter /), 
tsc{h) increases. This is related to the experimental observation [|l| that the irreversibility 
line shifts to higher temperatures with increasing dose of irradiation with heavy ions. 

As temperature drops below tsc{h), x remains at the value it had at the transition. 
There is now a macroscopic occupancy of the lowest energy state at = and k = 0. As is 
well known, condensation into this state is a meaningful concept only if the state is extended. 
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Condensation into a localized lowest energy state would imply a diverging Abrikosov ratio 
{(3a) below the transition, making the transition impossible. It is a special feature of this 
problem that the lowest lying state must indeed be extended for certain range of impurity 
concentrations. The density of states, Eq. 4, changes from being infinite at = when 
/ < 2, to being zero when / > 2. Thus, for fields and impurity concentrations such that 



parameter / < 2, true extended states (which always exist in the LLL [jT4[) must lay at the 
bottom of impurity band, since the number of states there diverges. The change of behavior 
in the density of states at / = 2 could be caused by the fact that the mobility edge shifts 
to positive energies at some /o > 2, leaving spread-out but localized states at \^ = 0, which 
now becomes the tail of the distribution. Numerical diagonalization studies indicate that 



mobility edge is indeed located near the band center for / > 4 [T^ . Thus, strictly speaking, 
our model is appropriate for / < /o- However, even for / above but close to /o, which is 
often the case for fields and concentrations of experimental interest, the states at = are 
still near mobility edge and will appear extended in a finite size sample. On this basis, we 
expect that useful information about the transition can still be obtained within our model. 

With these cautionary remarks in mind, the natural order parameter is the thermal 
average of the component of '?/'„(r) corresponding to the eigenvalue with V = Q and A; = 0. 
This is (V^co) = {NNi{ga\t=t^^[h) — So)/ {Pa)Y^'^- The disorder-averaged value of the field is 
{i^nif)) = 0, due to random phases of the state (t)v=oif)- Under the assumption that the 
lowest state is extended through the sample, 10^=0(^)1^ ~ the Edwards-Anderson order 



parameter |jT6[ qea = |('0n(r'))P then equals 



qEA = l-TT {9a\t=ts,ih) - Qa) (8) 

\Pa) 

below tsc{H), and is zero above. Thus, qea oc (tsc(^) — t)"^^, with the exponent P = 1/2. 
The free energy below tsc{h) is 

F gl , 1 r 



I l^^ POO 

+ 27-/^70 PA^)d^^''Mk)+9xV). (9) 



NNlT {13a) 

4 Correlation lengths and dynamical exponent 

To calculate the exponents that determine the divergence of correlation lengths parallel and 
perpendicular to the field we first note that from Eq. 6 and the definition of critical line it 
follows 



[g^ + X)[l + [ dkj 

4 J-i Jo 



(grieik) + gxV) {g^e{k) + gxV + ga + x) 

= ga- ga\t=tsc(h)- (10) 

The integral in the last equation diverges for / < 5/2 as {g^ + x)^^^^"^ when the transition 
line is approached from above, and it is finite for / > 5/2. Thus, we obtain {ga + x) (x 
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{t - tscY^^f-^/^^ for / < 5/2 and + x) oc (t - for / > 5/2. The same behavior follows 
if the transition line is approached along the line of constant temperature. This determines 
the value of the exponent for the correlation length parallel to the field oc (t — tsdh))'^''^ '■ 

z^ii = 1/(2/ - 3), (11) 

for / < 5/2 and the classical value v\\ = 1/2 for / > 5/2. The concentration corresponding 
to / = 5/2 determines the effective upper critical dimension in the problem. We now turn to 
correlation length perpendicular to the field, ^± oc [t — tsc{h)]~'^^ and study the susceptibility 
associated with Edwards- Anderson order parameter, 



XEA{r-r') = {i^*{r)ipn{r')){ipn{r)'ipn{r')). (12) 
After expanding the field operators in the eigenbasis of random potential we obtain 

XEAif- r') = I dV,dV2dhdk2F{f- P, V,, V2) ^^g^ 

y {gne{ki) + gxVi + go, + x){g^e{k2) + gxV2 + ga + x) 

where the function F is the two-particle spectral density [11 



Fir- r', V2) = J2 ^(Vi - - V,)0*(f)0,(r')0,(r-')0*(r') (14) 



and (f)i{f) are the eigenstates of the random potential. If we now introduce V = (Vi + V2)/2 
and uj = [Vi — V2)/2, for V close to mobility edge and small {q,uj), the Fourier transform of 
F has a diffusive form |[1^, [1^ 

Fi,i,V,.)^ 'fj^yjl'lf,, , (15) 

where D{q^/uj) is the generalized "diffusion constant". Assuming this form for F(g, V,a;) 
and rescaling everything by the appropriate power of temperature in Eq. 12, we find: 

i/x = (16) 

Surprisingly, the scaling turns out to be isotropic, in spite of the strong anisotropy in the 
model. 

The expression for two-particle spectral density yields immediately another important 
result. Since the energy and the momentum variables appear always in combination cu/g^, 
the dynamical exponent has the value: 

z = 2. (17) 

This may also be found by directly calculating the dc conductivity along the field using the 
time-dependent version of the model, which yields azz oc ^y, the same as Aslamazov-Larkin 
result in three dimensions. 
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The dependence of the correlation length exponents on magnetic field through parameter 
/ is the consequence of power-law behavior of the density of states close to the bottom of 
the band. A closer look at the the density of states Pf{V) reveals it to be roughly constant 
except in a narrow region, typically less than 1% of total bandwidth, around = 0, where 
it either diverges of vanishes |13[. One might expect that taking the effective potential V^(r^ 



to have a finite range would tend to wash away this fine feature of the density of states. 
This suggests that the experimentally relevant situation corresponds to / = 2 in our model, 
when the density of states is fiat down to the lowest energies. Based on this observation one 
would expect that the observable values of the exponents are = z/|| = 1. This expectation 
is met by the experimental results of Ref. 5: ux = 1.0 ± 0.1, i^y = 1.1 ± 0.1. The value of 
dynamical exponent is also in agreement with the measured z = 2.2 ± 0.2 ^j. 

5 Magnetization and specific heat 

Magnetization per unit volume is found by differentiating the averaged free energy with 
respect to the magnetic field: 



^ea + tttt]- (18) 



Below the transition line this coincides with the usual mean-field result. Above the transition 
line qea = and from Eq. 10 it follows that at constant temperature close to the transition 
{da + x) (X [h — hsc{t)]^^^-^~^^'^^ when / < 5/2 and {ga + x) oc [h — hsc{t)] otherwise. Thus the 
magnetic susceptibility is a smooth function of the field at the transition for / < 2, but has an 
upward cusp for 2 < / < 5/2 and a discontinuity for /> 5/2. The size of this discontinuity 
depends on the location of the transition in the H — T diagram. Differentiating the free 
energy twice with respect to temperature one obtains the specific heat. It is straightforward 
to show that at the transition it behaves the same way as susceptibility; smooth for / < 2, 
has a cusp for 2 < / < 5/2 and has the usual discontinuity for /> 5/2. More precisely, both 
magnetic susceptibility and specific heat behave as [t — tsc{h)]~" for 3/2 < / < 5/2, where 
a = (/ — 5/2)/(/ — 3/2). The behavior of the specific heat, order parameter and correlation 
length in our model is related to the one obtained from 0(2A^) vector model in the limit 
N —>■ oo and in the effective dimension -De// = 2/ — 1. 



6 Perturbations 

Having a solvable version of the general theory (1), we are in position to study small pertur- 
bations around it and check their relevance for the critical behavior. First, let us assume that 
the strength of disorder A is also allowed to fiuctuate from one columnar defect to another, 
i. e., the random potential in Eq. 1 is taken to be 

(19) 
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with the random variable Aj distributed according to the distribution: 



P{X^) = ^, (20) 

for |Aj — A| < w, < w < X, and zero otherwise. This random potential, being still spatially 
uncorrelated, allows an exact calculation of the density of states |T^. One can show however 



T3| , that the power-law behavior of the density of states at low energies is not changed by 
allowing finite width w, which implies the same critical exponents as already obtained. This 
perturbation is therefore irrelevant. 

A more complicated extension of the model is to assume a particular form of three- 
dimensional disorder represented by a random potential in Eq. 1: 

XnJ2^(^-^^)^ (21) 

i 

where the effective strength of disorder A„ is the same for all defects within the n-th layer, 
but fluctuates from one layer to another according to a Gaussian distribution 

P[{A„}] (X exp - J2{X^ - Xf/w\ (22) 

n 

where the width of distribution w is a small parameter, w << A. Since the model is not fully 
three, but rather 2-1-1 dimensional, one can still study the phase transition by the method 
described here in combination with the replica formalism |16|]. Without going into details. 



the result is that the critical behavior is not changed by a small w, and the only effect of 
additional disorder along the field is a decrease in transition temperarature. We conclude 
that this perturbation is also irrelevant for the critical behavior. 

As emphasized earlier, the replacement of disorder potential representing disorder by a 
sum of delta-functions is sensible only if the other two length scales in the problem, magnetic 
length and the average separation between the defects, are much larger than the average 
diameter of the columns. The latter condition obviously breaks down for higher doses of 
irradiation, when there is a considerable probability of overlap between the columnar defects. 
In fact, at higher densities of defects the trend of increase of the transition temperature with 
defect concentration should be reversed. To see this, consider first a small density of defects. 
Then our model is applicable, and as we already noted, at fixed magnetic field the transition 
temperature given by eq. 7 increases with increasing parameter /. This is because larger / 
shifts the density of states towards higher energies, decreasing the integral in denominator 
in the Eq. 7. For high defect densities the model ceases to be useful, but it is still possible to 
qualitatively understand what happens by considering the limit of very large concentrations. 
For very large density of defects, the effect of randomly distributed finite width scatterers is 
similar to having a potential which is almost constant in space. In such a potential the LLL 
is again nearly degenerate, only the average energy is shifted by some amount from where it 
was without any external potential. But this being the case, thermal fluctuations are again 
able to suppress the transition temperature because of the dimensional reduction! Thus, 



8 



for large densities of defects transition temperature must be decreasing with the density. 
This suggests that there is a field dependent optimal defect concentration at which the 
transition temperature reaches its maximum value, still below the mean-field Tc2{H). This 
would explain the observation of saturation of transition temperature with increase of dose 
of irradiation in ref. 16. 

7 Summary and conclusions 

We have studied the high-field superconducting glassy transition induced by columnar disor- 
der. It was argued that in certain range of magnetic fields and for not too large concentrations 
of defects the random potential that represents disorder may be assumed to be uncorrelated in 
space. By assuming the self-interaction of the fluctuating superconducting order-parameter 
to be infinitely ranged in directions perpendicular to the field we defined a model which 
describes the crossover regime towards the true glassy critical behavior, and which can be 
solved exactly. Within this model we have calculated the transition line, Edwards-Anderson 
order parameter, magnetization, specific heat, critical exponents for the correlation lengths 
and the dynamical exponent. The critical behavior of the model is shown to be isotropic, and 
the values of the critical exponents are in agreement with existing measurements. Several 
extensions of our model are considered, and some are shown to be irrelevant for the critical 
behavior. We presented a qualitative arguments for the behavior of transition temperature 
with density of defects, and argued for the existence of optimum concentration of defects at 
which the transition temperature reaches its maximum. 

The crucial simplifications of the general GL theory with columnar disorder in Eq. 
1 which facilitated the exact solution of the model were the assumption of delta-function 
scatterers for representation of disorder and the neglect of weak SCDW correlations, i. e., 
101^ interaction term was replaced by interaction infinitely long-ranged within layers. The 
first approximation is not a drastic one for small densities of defects; allowing for some range 
of random scatterers is likely to wash away the power-law behavior of density of states at 
low energies and make the critical behavior analogous to the / = 2 case for our model. This 
picture is in agreement with the measured values of critical indices for correlation lengths [Q. 
Neglect of SCDW correlations represents a more serious approximation, which necessarily 
breaks down sufficiently close to the critical temperature. The model solved here is thus 
best understood as a description of the crossover region before the true critical (glassy) 
behavior is reached. It is a disordered equivalent of spherical model for spin systems. Strong 
disorder however, guarantees a rather wide crossover region where our model is useful. This 
is supported by experimental observation of nearly isotropic scaling [^], while the general 
model in Eq. 1 would be most likely to exhibit an anisotropic critical exponents. Also, for 
calculations of non-universal quantities, such as the transition line in if — T phase diagram, 
the present model augmented with more realistic potential V{r) should suffice. It is therefore 
a useful step in quantitative understanding of superconducting transition in presence of 
correlated disorder in magnetic field. 
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